METHOD OF APPLYING PARAMETRIC OSCILLATORS TO MODEL DIELECTRIC 
FUNCTIONS 

This Application Claims benefit of Provisional Application Serial 
5 No, 60/471,753 Filed 05/20/03. 

TECHNICAL FIELD 

The present invention comprises methodology for modeling 
10 data, and more particularly a method for segmentally modeling 
Real and Imaginary Parts of Optical Constants or Dielectric 
Functions using Kramers-Kronig (K-K) Consistant Oscillators that 
have substantially definite start and end points. 

15 BACKGROUND 

Empirical investigation of substrate systems allows 
determination of the thicknesses of/ and associated optical 
property component values, (eg. refractive index (n) and 
20 extinction coefficient (k), or mat-hkematically equivalent real 

(el) and imaginary (e2) complex dielectric function values), for 
films thereon. 

Optical property component values for substrate systems can 
25 be empirically obtained by many methods, such as those utilizing 
substrate system film effected changes caused in reflected and/or 
transmitted light beam intensities, and those utilizing substrate 
system film effected changes in ellipsometric PSI and DELTA 
values. However obtained, resulting optical property component 
30 values, when plotted as a function of the wavelength, (ie. photon 
energy, which comprises of a beam of electromagnetic radiation 
Utilized in the Investigation), typically present as a rather 
irregular plot, which irregular plot can not be accurately 
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represented by a simple analytical mathematical function over the 
range of wavelengths, (ie, photon energies), typically of 
interest, (eg, zero (0.0) ev to six (6.0) ev, and above). 

5 To understand the present invention, it is helpful to 

realize that the index of refraction of a substrate system film 
represents the ratio of the speed electromagnetic radiation 
traveling in vacuum to that in said film material, and that it is 
typically a complex number, (N = n + ik), where N is the Index of 

10 Refraction, n is termed the refractive index and k is termed the 
extinction coefficient. An alternative representation can be 
arrived at by the relation (e = N'^2), and is termed the 
Dielectric Function, (ie. e = el +ie2), where (el = n"2 - k''2) 
and (e2 = 2 x n x k). It is also noted that the photon energy 

15 of a beam of electromagnetic radiation is related to the 

frequency thereof by the equation (E = h x F), where E is energy, 
(typically in electron volts ev), h is Plank's Constant, and F is 
Frequency. As well. Frequency, in electromagnetic radiation, is 
related to Wavelength by the relationship (C = F x W), where C is 

20 the speed of light and W is wavelength, hence (E = h x C/W) . In 
addition, electrical properties of a material are related to the 
optical properties by known relationships. For instance, 
electrical conductivity is found to be proportional to e2. 

25 One approach to modeling the Imaginary Part (e2) of a 

Dielectric Function is to identify some number of points 
thereupon and fit a polynomial thereto. Applying a Kramer-Kronig 
(K-K) Integration to the Imaginary Part (e2) provides the Real 
Part (el). This is discussed in an Article titled "High 

30 Precision UV-Visible-Near-IR Stokes Vector Spectroscopy", Zettler 
et al.. Thin Solid Films, 234 (1993). It is also important to 
understand that the real and imaginary components of the 
Refractive Index (ie. n and k), are not independent of one 
another. As principal of causality applies, known as the 
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Kramer-Kronig, (K-K), relationship, and theoretically, allows 
determining k(E} when n(E) is known, and vice-versa. 

Continuing, other Dielectric Function Models have been 
5 developed for semiconductors. Adachi has developed a (K-K) 

consistent model suitable for below energy band-gap calculations 
that has been used to describe ternary alloys by interpolating 
parameters between binary endpoints, (see S. Adachi, "Optical 
dispersion relations for 6aP, GaAs, GaSb, InAs InSb, 

10 Al(x)Ga(l-x)As, and In(l-x)As(y)P(l-y) , Appl. Phys., Vol 66, p. 

6030, 1989. The ability to interpolate parameters produces much 
more physically realistic dielectric models than simply averaging 
dielectric functions of the end-point binarjr materials. For 
terniary alloy models, the strong Critical Point (CP) structure 

15 of semiconductors causes simple averaging schemes to produce 

doubled (CP) structures at the energies of binary endpoints, (see 
P.O. Snyder, J. A. Woollam, S.A. Alterovitz and B. Johs, "Modeling 
Al(x)Ga(l-x)As Optical Constants as Functions of Composition", J. 
Appl. Phys., Vol. 68, p. 5925, 1990). Forouhi and Bloomer have 

20 provided a dielectric model for semiconductors which is (K-K) 
consistent and has a very small set of coefficients, (see A.R. 
Forouhi and I. Bloomer, "Optical Properties of Crystaline 
Semiconductors and Dielectrics", Phys. Rev. B, vol. 38, p. 1865, 
1988 and Patent No. 4,905,170). The Forouhi and Bloomer model, 

25 however, has been found to have insufficient flexibility to fit 

existing dielectric functions accurately enough for ellipsometric 
modeling. Oscillator ensembles, (eg. harmonic and Lorentz), 
have been used to describe the above energy-gap behavior of some 
semiconductors and the AlGaAs alloy system, (see a paper by M. 

30 Erman, J.B. Theeten, P.Chambon, S.M. Kelso and D.E. Aspnes, 

titled "Optical Properties and Damage Analysis of GaAs Single 
Crystals partly Amorphised by ion Implantation", J. Appl. Phys. 
vol. 56, p. 2664, 1984; and a paper by H.D. Yao, P.G. Snyder and 
J .A. Woollam, titled "Temperature Dependence of Optical 
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Properties of GaAs", J. Appl. Phys., vol. 70, p. 3261, 1991; and 
a paper by F. Terry Jr., "A Modified Harmonic Oscillator 
Approximation Scheme for the Dielectric Constants of 
Al(x)Ga(l-x)As", J. Appl., vol 70, p. 409, 1991). These models 
5 have been used to fit measured ellipsometric data, however, they 
are incapable of describing direct *energy-band-gap spectral 
regions, and they require extra fictitious oscillators to fill in 
the absorbtion between (CP's). Kim and Garland et al. have 
developed a (K-K) consistent model that can adequately describe a 

10 semiconductor dielectric function above, below and through the 

fundamental direct energy-gap, and this model has been applied to 
the AlGaAs alloy system, (see C.C. Kim, J.W. Garland, H. Abad and 
P.M. Raccah, **Modeling the Optical Dielectric Function of 
Semiconductors: Extension of the Critical- Point-Parabolic-Band 

15 Approximation, Phys Rev. B, vol 45, p. 11749, 1982; C.C. Kim, 
J.W. Garland, H. Abad and P.M. Raccah, "Modeling the Optical 
Dielectric Function of the Alloy System Al(x)Ga(l-x)As*', Phys. 
Rev. B, vol. 47, p. 1876, 1993). This model can accurately 
describe the dielectric function and higher order derivatives. 

20 However, to determine required internal coefficients a two stage 
fitting process is used. First (CP) energies and broadening are 
determined by fitting derivatives of the dielectric function, and 
then the remaining internal coefficients are determined with the 
energies and broadenings fixed. It will be appreciated that this 

25 model then requires that the dielectric functions exist before 
the model can be fitted. Furthermore, attempts at fitting all 
coefficients simultaneously (as necessary for direct ellipsometer 
data fitting), are unlikely to succeed because of the highly 
correlated nature of the functions internal to the model. Over 

30 part of the spectral range, the modeled imaginary part of the 
dielectric function results from the difference of internal 
function values one-hundred (ICQ) times larger than the final 
value. The internal coefficients are delicately balanced to 
produce the proper output, and small changes in (CP) energies can 
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cause large deviations in the model output. In addition, 
Lorentzian Broadening, which is known to be wrong for elements 
and compounds, is utilized in by Kim & Garland et al. work* 
While the Kim & Garland et al. work is very interesting, there 
5 remained need for an improved Parametric Model such a that 

taught in Patent No. 5,796,983, which is discussed in detail in 
the Disclosure and Detailed Description Sections of this 
Specification, along with other applicable Oscillator Structures, 
(eg. Gaussian, Marrow Lorentzian, Harmonic, TOLO, lonicl, 

10 Ionic2). As demonstrated in the Detailed Description Section of 
this Specification, typical conventional practice is to place 
Oscillator Structures in a Plot of the Imaginary Part of the 
Dielectric Function (e2) such that at each wavelength the sum of 
the contributions results in said Imaginary Part. Typically this 

15 requires placing Oscillator Structures under peaks, and 
« elsewhere, which placement is an acquired art. 

Need remains for a simple mechanical approach to placing 
Oscillator Structures in a Plot of the Imaginary Part of the 
20 Dielectric Function (e2) such that summation of contributions 

thereof at each wavelength results in said Imaginary Part of the 
Dielectric Function (e2). Said approach should utilize (K-K) 
consistant Oscillator Structures so that the Real part of the 
Dielectric Function results from an integration procedure. 
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niSCLOSURE OF THR INVEMTION 

TO begin, it is to be generally understood that the practice 
o£ ellipsometry involves causing polarized electromagnetic 
5 radiation to interact with a sample, and determining sample 

associated ellipsometr ic PSI (Y) and DELTA (A) values, where 
PSI (y) and DELTA (A) a" related by: 
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p = rp/rs = Tan(T) exp(iA) 



Ellipsometric PSI and DELTA correspond to change in a ratio of 
orthogonal beam components (typically (rp/rs) where "p" and "s" 

15 components are in a "p" plane defined by the beam and a 
perpendicular to the sample surface, and in an "s" plane 
perpedicular to said "p" plane and parallel to the sample 
surface), and change in phase between the "p" and "s" components 
respectively. Dielectric Functions for a sample, which consist 

20 of real (el) and imaginary (e2) components, can be derived from 
ellipsometric PSI and DELTA values via a regression procedure. 

It is noted that it is prefered practice to fit oscillator 
structures to the Imaginary Part (e2) of the Dielectric Function 

25 such that summation of contributions therefrom at each wavelength 
provide a fit. Typical practice is to position oscillator 
Structure peaks under peaks in the (e2) plot. The disclosed 
invention breaks with that convention and simply places a 
plurality or multiplicity of typically, though not necessarily, 

30 similar oscillator structures at identified wavelength locations 
which do not necessarily have any definite relationship to the 
shape of the Dielectric Function. Some will be positioned so 
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that their peaks are located near (e2) plot peaks, but most are 
not so positioned. The spacing of the oscillator structures 
under present invention teachings can be equal, or can be 
substantially random, or can involve more oscillator structures 
being placed in wavelength regions near peaks in the (e2) plot. 
Any functional arrangement is within the scope of the present 
invention, but the goal is not to model a plot with a minimum 
number of oscillator structures possible as is the goal in 
conventional procedures. Instead the goal is to make the 
procedure mechanical so that essentially anyone, regardless of 
training and experience, can practice it successfully after a 
very short period of instruction. 

A disclosed invention method of segmentally modeling real 
and imaginary parts of dielectric functions with Kramers-Kronig 
(K-K) consistency, then comprises the steps of: 

a) providing an imaginary part of a dielectric function and 
dividing it into a plurality of segements; 

b) fitting each segement with an appropriate approximating K-K 
consistent mathematical construction, the mathematical 
constructions being continuous at intersection between adjacent 
segements, each mathematical construction comprising a selection 
from the group consisting of: 

a mathematical equation; and 

a summation of oscillator structures. 

Where a mathematical equation is fit to the imaginary imaginary 
part of a dielectric function it can be of a polynomial form of 
degree (n-1), where (n) is sum of the number of intersections 
between segements plus a beginning and an end node. It is 
disclosed, however, that the technique of fitting the imaginary 



part of a Dielectric Function vith a polynomial of degree one 
less than the number of points being fit, (eg* a Cubic Spline as 
discussed in "Applied Numerical Methods for Digital Computation", 
James et al.. Harper Row (1985)), is known, followed by applying 
a K-K Integration to provide the real part of the Dielectric 
Function. In contrast, the disclosed invention provides a 
separate mathematical construction for each segement between 
identified points. 

A more relevant presently disclosed method of segmentally 
modeling real and imaginary parts of dielectric functions with 
Kramers-Kronig (K-K) consistency, then comprising the steps of: 

a) providing an imaginary part of a dielectric function over 
spectroscopic range, and dividing said spectroscopic range into a 
plurality of segements; 

b) fitting each said segement in said spectroscopic range with 
an approximating K-K consistent oscillator structure, said 
approximating oscillator structure in each said segement 
beginning and ending at the start and end of said respective 
segement, such that a summation of contributions from said 
oscillator structures present at each point within said 
spectroscopic range approximates said imaginary part of said 
dielectric function, and via Kramers-Krdnig (K-K) consistency, 
also the real part of said dielectric function 

A specific embodiment of said more relevant method provides that 
the segements are of equal spectroscopic range lengths and at 
least one of said K-K consistent oscillator structures is 
preferably triangular shaped; the start and end of all oscillator 
structures, except the start of the first and end of the last, 
being positioned at the same spectroscopic point as are peaks of 
immediately adjacent oscillator structures. 
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Another specific embodiment of said more relevant method provides 
that at least one segement in said spectroscopic range is of a 
different length than other segements in said spectroscopic range 
and at least one of said K-K consistent oscillator structures is 
5 preferably triangular shaped; the start and end of all oscillator 
structures, except the start of the first and end of the last, 
being positioned at the same spectroscopic points as are peaks of 
immediately adjacent oscillator structures « 

10 A more detailed disclosure of a present invention method of 

seqmentally modeling real and imaginary parts of dielectric 
functions using Kramers-Kronig (K-K) consistent oscillators, 
comprises the steps of: 

15 practicing steps a and b in either order; 

a) providing experimentally obtained data for real and 
imaginary parts of a dielectric function vs. wavelength for a 
sample comprising a transparent thin film on a substrate, over a 

20 determined range of wavelengths; 

b) providing a mathematical model of said sample which 
comprises a parameter coresponding to the thickness of the 
transparent thin film and comprises parameters corresponding to a 

25 pole amplitude and location at a wavelength beyond the lower 
wavelength, high energy, extent of said determined range of 
wavelengths, then by a global fitting procedure evaluating 
parameters, including the thickness of the transparent thin film 
and the location and amplitude of the pole in said mathematical 

30 model, utilizing the data corresponding to the real part of the 
dielectric function; 

with transparent thin film thickness evaluated in step b. 
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proceeding to practice steps c, d and e sequentially: 

c) defining a wavelength range segement length such that the 
sum of (n + 0.5) wavelength range segements exactly span the 

5 determined wavelength range, and beginning at one end of said 

determined wavelength range placing a (K-K) consistent oscillator 
which comprises an amplitude parameter and begins and ends at the 
extents of the first segement with its peak midway therebetween, 
then performing a point by point fit to the Imaginary part of the 
10 dielectric function data over wavelengths in said first 

wavelength range segement such that (K-K) consistent oscillator 
defining parameters are evaluated; 

d) placing a second (K-K) consistent oscillator which begins at 
15 a wavelength at which the first (K-K) consistent oscillator peaks 

and ends one wavelength range segement length therefrom and has a 
peak midway therebetween, then performing a point by point fit to 
the imaginary part of the dielectric function data over said 
first and second wavelength range segements such that (K-K) 
20 consistent oscillator defining parameters in said first (K-K) 
consistent oscillator are re-evaluated and oscillator defining 
parameters in said second (K-K) consistent oscillator are 
evaluated; 

25 each of the remaining (n - 2) wavelength range segements, 

sequentially; 

e) placing a (K-K) consistent oscillator which begins at a 
wavelength at which the Just prior (K-K) consistent oscillator 

30 peaks and ends one wavelength range segement length therefrom and 
having a peak midway therebetween, then performing a point by 
point fit to the imaginary part of the dielectric function data 
over said wavelength range segements which are fitted with (K-K) 
consistent oscillators, such that oscillator defining parameters 
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in previously evaluated (K-K) consistent oscillators are 
re-evaluated and oscillator defining parameters in the added 
oscillator are evaluated; 

such that at each wavelength over the determined range of 
wavelengths the sum of the contributions of each evaluated (K-K) 
consistent oscillator approximates the magnitude of the imaginary 
part of the dielectric function. 

The just recited method assumes use of equal length 
wavelength range segements, but said method of segementally 
modeling real and imaginary parts of dielectric functions with 
Kramers-Kronig (K-K) consistent oscillators could involve 
utilizing^a^^^^^ length wavelength range segements, and 

comprise the steps of: 

practicing steps a and b in either order; 

a) providing experimentally obtained data for real and 
imaginary parts of a dielectric function vs. wavelength for a 
sample comprising a transparent thin film on a substrate, over a 
determined range of wavelengths; 

b) providing a mathematical model of said sample which 
comprises a parameter coresponding to the thickness of the 
transparent thin film and comprises parameters corresponding to a 
pole amplitude and location at a wavelength beyond the lower 
wavelength, high energy, extent of said determined range of 
wavelengths, then by a global fitting procedure evaluating 
parameters, including the thickness of the transparent thin film 
and the location and amplitude of the pole in said mathematical 
model, utilizing the data corresponding to the real part of the 
dielectric function; 
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with transparent thin film thickness evaluated in step b, 
proceeding to practice steps c, d and e sequentially: 

c) defining n wavelength range segement lengths such that the 
sum of said n wavelength range segements plus half the length of 
the last wavelength range segement at one end of saiid wavelength 
range exactly spans the determined wavelength range, and 
beginning at the opposite end of said determined wavelength range 
placing a (K-K) consistent oscillator which comprises an 
amplitude parameter and begins and ends at the extents of the 
first segement with its peak midway therebetween, then performing 
a point by point fit to the imaginary part of the dielectric 
function data over wavelengths in said first wavelength range 
segement such that (K-K) consistent oscillator defining 
parameters are evaluated; 

d) placing a second (K-K) consistent oscillator which begins at 
a wavelength at which the first (K-K) consistent oscillator peaks 
and ends the second wavelength range segement length therefrom 
and has a peak midway therebetween, then performing a point by 
point fit to the imaginary part of the dielectric function data 
over said first and second wavelength range segements such that 
(K-K) consistent oscillator defining parameters in said first 
(K-K) consistent oscillator are re-evaluated and oscillator 
defining parameters in said second (K-K) consistent oscillator 
are evaluated; 

for each of the remaining (n - 2) wavelength range segements, 
sequentially; 

e) placing a (K-K) consistent oscillator which begins at a 
wavelength at which the just prior (K-K) consistent oscillator 
peaks and ends at the nth wavelength range segement length 
therefrom and having a peak midway therebetween, then performing 
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a point by point fit to the imaginary part of the dielectric 
function data over said wavelength range segements which are 
fitted with (K-K) consistent oscillators^ such that oscillator 
defining parameters in previously evaluated (K-K) consistent 
oscillators are re-evaluated and oscillator defining parameters 
in the added oscillator are evaluated; 

such that at each wavelength over the determined range of 
wavelengths the sum of the contributions of each evaluated (K-K) 
consistent oscillator approximates the magnitude of the imaginary 
part of the dielectric function. 

In either of the embodiments recited above, said method can 
further comprise step f: 

f ) using starting values of parameters previously evaluated 
performing a global fitting procedure onto both real and 
imaginary parts of the dielectric function to re-evaluate 
parameters, including the thickness of the transparent thin film, 
the location and amplitude of the pole in said mathematical 
model, and the "N" (K-K) consistent oscillator defining 
parameters utilizing the data corresponding to the real and 
imaginary parts of the dielectric function, with a constraining 
limitation that said pole location is required to remain located 
at a wavelength beyond the lower wavelength, high energy, extent 
of said determined range of wavelengths. 

It should be appreciated that in either version, the pole 
location and its amplitude at a wavelength beyond the lower 
wavelength, high energy, extent of said determined range of 
wavelengths, can be fixed during steps c, d and e, or the pole 
location and/or its amplitude can be re-evaluated along with 
previously evaluated oscillator defining parameters in previously 
evaluated (K-K) consistent oscillators during steps c, d and e. 
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but wherein said pole location is required to remain located at a 
wavelength beyond the lower wavelength, high energy, extent of 
said determined range of wavelengths. The later constraint is 
necessary to prevent the Pole form entering large influence 
within the determined range of wavelengths. This is to be 
appreciated in the light that the purpose of said Pole is to 
provide a catch-all means to account for all influences beyond 
the determined range of wavelengths. 

It is to be appreciated that the (K-K) consistent 
oscillators can be of a preferred tiangular shape, or each (K--K) 
consistent oscillator can be Gaussian or of a narrow Lorentzian 
shape, or constructed from at least one polynomial on at least 
one side of the peak value thereof. 

In addition, said method, for each of the "n" the wavelength 
range segements, can involve letting the segement length thereof 
float and be fit along with other (K-K) consistent oscillator 
parameters . 

Also, the global and point by point fits are based on a 
square error minimization criteria. 

A modified method of segementally modeling real and 
imaginary parts of dielectric functions with Kramers-Kronig (K-K) 
consistent oscillators, comprises the steps of: 

practicing steps a and b in either order; 

a) providing experimentally obtained data for real and 
imaginary parts of a dielectric function vs. wavelength for a 
sample comprising a transparent thin film on a substrate, over a 
determined range of wavelengths; 
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b) providing a mathematical model of said sample which 
comprises a parameter coresponding to the thickness of the 
transparent thin film and comprises parameters corresponding to a 
pole amplitude and location at a wavelength beyond the lower 
wavelength, high energy, extent of said determined range of 
wavelengths, then by a fitting procedure evaluating parameters, 
including the thickness of the transparent thin film and the 
location and amplitude of the pole in said mathematical model, 
utilizing the data corresponding to the real part of the 
dielectric function. 

With transparent thin film thickness evaluated in step b, said 
modified method then proceeds to with steps c, d and e 
sequentially: 

c) defining a wavelength range segement length such that the 
sum of (n) wavelength range segements exactly span the determined 
wavelength range, and beginning centrally in said determined 
wavelength range placing a first (K-K) consistent oscillator 
which comprises an amplitude parameter and begins and ends at the 
extents of the first segement with its peak midway therebetween; 

for each of the remaining (n - 1) wavelength range segements, on 
ether side of the central peak of the first (K-K) consistent 
oscillator; 

d) placing a (K-K) consistent oscillator which begins at a 
wavelength at which the just centrally prior (K-K) consistent 
oscillator peaks and ends one wavelength range segement length 
therefrom and having a peak midway therebetween; 

e) performing a fit to the imaginary part of the dielectric 
function data over said wavelength range segements which are 
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fitted with (K-K) consistent oscillators, such that oscillator 
defining parameters in the (K-K) consistent oscillators are 
evaluated; 

such that at each wavelength over the determined range of 
wavelengths the sum of the contributions of each evaluated (K-K) 
consistent oscillator approximates the magnitude of the Imaginary 
part of the dielectric function. 

The just recited modified method assumes that the wavelength 
range segement lengths are equal. Another modified method of 
segementally modeling real and Imaginary parts of dielectric 
functions with Kramers-Kronlg (K-K) consistent oscillators allows 
for the wavelength range segements to vary, and comprises the 
steps of: 

practicing steps a and b In either order; 

a) providing experimentally obtained data for real and 
Imaginary parts of a dielectric function vs. wavelength for a 
sample comprising a transparent thin film on a substrate, over a 
determined range of wavelengths; 

b) providing a mathematical model of said sample which 
comprises a parameter corespondlng to the thickness of the 
transparent thin film and comprises parameters corresponding to a 
pole amplitude and location at a wavelength beyond the lower 
wavelength, high energy, extent of said determined range of 
wavelengths, then by a fitting procedure evaluating parameters, 
including the thickness of the transparent thin film and the 
location and amplitude of the pole In said mathematical model, 
utilizing the data corresponding to the real part of the 
dielectric function; 
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with transparent thin film thickness evaluated in step 
proceeding to practice steps c, d and e sequentially: 



c) defining (n) wavelength range segement lengths such that the 
5 sum of said (n) wavelength range segements exactly span the 

determined wavelength range^ and beginning centrally in said 
determined wavelength range placing a first (K-K) consistent 
oscillator which comprises an amplitude parameter and begins and 
ends at the extents of the first segement with its peak midway 
10 therebetween; 

for each of the remaining (n - 1) wavelength range segements, on 
ether side of the central peak of the first (K-K) consistent 
oscillator; 

15 

d) placing a (K-k) consistent oscillator which begins at a 
wavelength at which the just centrally prior (K-K) consistent 
oscillator peaks and ends one wavelength range segement length 
therefrom and having a peak midway therebetween; 

20 

e) performing a. fit to the imaginary part of the dielectric 
function data over said wavelength range segements which are 
fitted with (K-K) consistent oscillators, such that oscillator 
defining parameters in the (K-K) consistent oscillators are 

25 evaluated; 

such that at each wavelength over the determined range of 
wavelengths the sum of the contributions of each evaluated (K-K) 
consistent oscillator approximates the magnitude of the Imaginary 
30 part of the dielectric function. 

It is noted that said alternative procedures could involve 
selecting the wavelength range segement length to be the entire 
wavelength range, and involve use of only one centrally placed 
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oscillator. Further, it could involve use of three, or five etc. 
odd number oscillators. If fact, a procedure could be to begin 
with say one oscillator and perform a fit, and see if the fit is 
good enough. If not, additional wavelength range segements and 
oscillators can be added until a fit results in acceptable 
results. 

The pole location and its amplitude at a wavelength beyond 
the lower wavelength, high energy, extent of said determined 
range of wavelengths can be fixed during step e, or the pole 
location and/or its amplitude are re^evaluated along with 
previously evaluated oscillator defining parameters in the 
(K-K) consistent oscillators during step e, but wherein said pole 
location is required to remain located at a wavelength beyond the 
lower wavelength, high energy, extent of said determined range of 
wavelengths. 

The (K-K) consistent oscillators can be of tiangular shape, 
be constructed from at least one polynomial on at least one side 
of the peak value thereof, be of Gaussian, Narrow Lorentzian, 
Harmonic, TOLO, lonicl, Ionic2 shape, or be of any functional 
shape or construction. (See the Detailed Description and 
Drawings of this Specification, and the J. A. Woollam WVASE Manual 
which is incorporated by reference herein, for insight to the 
various identified Oscillator structures). 

It is again noted that the fits are preferably based on a 
square error minimization criteria such as, but not limited to, 
the Levenberg-Marquard approach. 

MATERIAL PREVIOUSLY DISCLOSED IN PATENT NO. 5^796.983 
INCLUDED TO PROVIDE INSIGHT TO MATHEMATICAL CONSTRUCTIONS 
COMPRISING NON-CONVENTIONAL OSCILLATOR STRUCTURES 
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Patent No. 5,796,983 discloses that problem with prior 
modeling approaches which utilize the summation of the effects of 
a plurality of Oscillator Structures to provide approximate 
values for (el) and/or (e2) as a function of photon energy, is 
that conventional Oscillator Structures do not provide accurate 
values of (el) and/or (e2) at all values of photon energy. While 
a Lorentzian or Gaussian Oscillator Structure might be capable of 
providing a good (MSB) fit to a "Peak" region in a plot of (e2) 
with respect to photon energy, for instance, unless a "Tail" 
region projection from said "Peak" region in said plot is 
precisely Mathematically Lorentzian or Gaussian, data points 
removed from the region of the Peak will not be well approximated 
thereby. That is, the shape of the "Tail" regions of the 
Lorentzian or Gaussian Oscillator Structure is dictated by the 
shape of the peak region because of a limited number of degrees 
of freedom available. In some cases the sum of a plurality of 
Lorentzian or Gaussian Oscillator Structures centered at various 
central energies, with various amplitudes and broadenings can be 
utilized to provide improved fit over a range of photon energies. 
However, even where this approach is utilized, the lack of a 
sufficient number of degrees of freedom in an Oscillator 
Structure based Model to allow shaping the Peak region of an 
Oscillator Structure independently of the shape of the Tail 
regions, normally prevents achieving a good (MSB) fit over a 
relatively large range of photon energies. 

To overcome the identified shortcommings in prior Oscillator 
Structure based Parametric Models for use in calculating 
Dielectric Function Component values, given a photon energy, and 
to meet the First Requirement recited in the Background Section, 
(ie. that, a Dielectric Parametric Function Model provide 
sufficient flexibility to fit essentially any empirically derived 
data set). That is. Oscillator Structures should preferably 
allow sufficient degrees of freedom to allow "Peak" and "Tail" 
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regions thereof to be mathematically modeled essentially 
independently. The preferred embodiment of such an Oscillator 
Structure provides that Finite Order Polynomials, (providing 
terms of any Order), and possibly Finite Magnitude 
5 Discontinuities, should be combined to form "Oscillator 
Structures". A typical 983 Patent Oscillator Structure is 
comprised of two Finite Order Polynomials (Fl and F2), to the 
left of a Center Point region, said two Finite Order Polynomials, 
(Fl and F2), being sequentially continuous in zero (0), first (1) 

10 and preferably second (2) derivatives at their point of merger; 
said two polynomials, (Fl and F2), being possibly followed by a 
Finite Magnitude Discontinuity; with the resulting foregoing 
described partial Oscillator Structure being followed by two 
additional Finite Order Polynomials, (F3 and F4), to the right of 

15 the Finite Magnitude Discontinuity, which Finite Order 

Polynomials, (F3 and F4), are sequentially continuous in zero 
(0), first (1) and preferably second (2) derivatives at their 
point of merger. It is also mentioned that a variation of a 
basic Oscillator Structure provides that some of the Finite Order 

20 Polynomials can be eliminated or made to contain only various 
Order Terms, (eg. that is contain any of zero, first, second, 
third etc. Order terms, and possibly not others), and still 
provide a functional Oscillator Structure. For instance, there 
might be only one Finite Order Polynomial on one or both sides of 

25 an Oscillator Structure Center Point, and a Finite Magnitude 

Discontinuity might, or might not, be present and might, or might 
not, be preceded or followed by a Finite Order Polynomial, for 
instance, in any case, each oscillator structure, (assuming more 
than one such Oscillator Structure is utilized), begins and ends 

30 at dependent variable (el) and/or (e2) values of zero (0), at an 
independent variable photon energy which is that corresponding to 
the Center Point of an Oscillator Structure, which can be that of 
the same, or another. Oscillator Structure. (Note if both the 
beginning and end of an Oscillator Structure are present at the 
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Center Point of the same Oscillator Structure, a Zero-Width 
Finite Magnitude Oscillator Structure "Pole" results, which is 
mathematically equivalent to a Finite Magnitude Discontinuity 
which begins at zero of a dependent variable scale at the Center 
5 Energy of the independent variable. This equivalency allows 
various approches to computer programming implementation). 

It is noted that Finite Order Polynomials and Finite 
Magntude Discontinuities are selected as components of Oscillator 

10 Structures because said mathematical entities can be subjected to 
a Convolution Integration, (a procedure involving, in the present 
case, the product of such a Finite Order Polynomial, with or 
without an associated Finite Magnitude Discontinuity; and a 
Broadening Factor, (preferrably Gaussian), which involves Error 

15 Functions and/or Exponentials, see supra). If only Finite Order 
Polynomials and Finite Magnitude Discontinuities are utilized to 
construct non-Zero-Width Oscillator Structures, analytic 
functions and published Integral Tables, (see the "Handbook of 
Mathematical Functions", Abromowitz and Stegun, Dover 

20 Publications, 1972 and "Table of Integrals, Series and Products, 
Gradshteyn and Ryzhik, Academic Press, 1980, which references are 
incorporated by reference herein), are available which allow 
relatively easy evaluation of required Oscillator Structure Model 
Convolution Integrals which arise from the application of an 

25 Oscillator Structure "Broadening Factor". 

As regards application in the presently disclosed invention, 
it is important to realize that a 983 Patent Oscillator Structure 
has definite beginning and ending points, and can be configured 
30 with single straight lines on either side of a peak, to form 
triangular shapes. 

Conventional practice typically involves inspection of 
Dielectric Function (e2) vs. photon energy to determine the 
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location(s) at which "Dependent Variable Peaks" exist. Use of 
plotted Dependent vs- Independent Variable data is useful in 
performing this inspection. The independent variable "Center 
Point" of an Oscillator Structure is then located at the photon 
energy corresponding to each Peak, with Polynomial Modeled Tails 
projecting to one or both sides thereof. An Oscillator Structure 
Finite Magnitude Discontinuity might also be positioned at a 
Center Point. (In the case where a Dependent Variable Peak 
corresponds to the Band-Gap Edge in a Semiconductor, Dielectric 
Function Component Magnitudes often demonstrate essentially 
Discontinuous "Jumps" over very small changes in Energy, at a 
photon energy corresponding to the edge of said Band-Gap. The 
inclusion of Finite Magnitude Discontinuities in the Model speaks 
to Sixth (6) Requirement recited in the Background Section of 
this Disclosure in which it is stated that a Model might provide 
physically relevant parameters therein, eg. a Bandgap absorption 
edge). Also, each Finite Width Oscillator Structure so 
positioned will begin and end at an independent variable (e2) 
value of zero (0), which zero (0) values are, (assuming more than 
one Oscillator Structure is present), located at the Center Point 
Energy of a positioned Oscillator Structure, which, as alluded to 
above, can be the same or another Oscillator Structure. This 
approach to locating the beginnings and ends of Oscillator 
Structures positioned on a plot of (e2) vs. photon energy, at the 
Center Point Energy of other Oscillator Structures, allows moving 
a Center Point Energy position in relation to other Center Point 
Energy positions without causing drastic effects in Model 
produced results at distal photon energies. This approach to 
locating lower and upper energy end-points of Oscillator 
Structures is absolutely essential to application of the present 
invention. 

With an "appropriately" estimated number of appropriately 
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constructed Oscillator Structure (s) positioned along a photon 
energy axis o£ an (e2) vs. photon energy plot^ a Finite Order 
Polynomial Coefficient and Finite Magnitude Discontinuity 
Magnitude defining Oscillator Structure Mathematical Model 
5 Parameter evaluation routine^ including applied Broadening 

Factors, can be applied to simultaneously provide values for all 
Oscillator Structure Mathematical Model Coefficients, based upon, 
for instance, a Mean-Square-Error (MSB) criteria. (A typical 
Coefficient evaluation routine might utilize a Levenberg-Marquard 

10 approach). Use of said Finite Order Polynomial Coefficients and 
Finite Discontinuity Magnitudes in constructing Oscillator 
Structures then allows evaluation of Oscillator Structure 
Mathematical Model Coefficients, in view of Convolution 
Integration applied Broadening Factors, which Mathematical Model 

15 Coefficients are utilized in evaluation of dependent variable, 

such as (e2) values, provided an independent variable value, such 
as photon energy. It will be appreciated that Oscillator 
Structure Mathematical Model Coefficient values for each 
Oscillator Structure utilized are all that have to then be 

20 available to allow calculation of a value of, for instance, (e2), 
given a value of, for instance, photon energy, when required • 
If, for example, the number of Oscillator Structures utilized is 
seven (7), and twelve (12) variables are found necessary to 
define each, then only eighty-four (84) Coefficient values must 

25 be stored. Thereafter, when a value of dependent variable (e2) 
for a given independent variable photon Energy level (ES) is 
desired, the various Finite Order Polynomials which are present 
in a range of photon energies in which the specific photon energy 
(ES) of interest is present are produced and evaluated at said 

30 photon energy (ES) value of interest, and the contributions of 
each Finite Order Polynomial are added together by a summation 
procedure. If necessary the value of a Finite Magnitude 
Discontinuity, (be it additive or subtractive ) , is included in 
the summation. 
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Nov, in use, one can try a number of "Mathematical Models". 
That is, it might be found that while seven (7) Oscillator 
Structures, each providing four (4) Polynomials, each of second 
(2) order, provide a good (MSB) fit over the range of Energies 
involved, some of the attendant Coefficients of the Mathematical 
Model for oscillator Structures are strongly "Correlated", such 
that change in one thereof is attended by a directly off-setting 
change in one or more of the others. That would be indicative 
of (e2) vs. photon energy Data which are not "strong" enough to 
justify the number of Coefficients requiring evaluation assumed 
into the Mathematical Model. One might then try eliminating one 
or more present Oscillator Structures or possibly one or more 
Finite Order Polynomials from one or more of the Oscillator 
Structures, or perhaps one might try reducing the number of 
Polynomial Order Terms in one or more Finite Order Polynomials 
present, and rerunning the Mathematical Model Coefficient 
evaluating program. If an equally good (MSB) fit is obtained, 
and less Mathematical Model Coefficient Correlation is present, 
the "Modified Mathematical Model" would be considered superior in 
a conventional sense. In use it has been found that far less 
than the above mentioned eighty-four (84) parameters are required 
to adequately characterize, for instance, the (e2) vs. photon 
energy Function of a Gallium-Arsonide Sample System. (Note, a 
highly corelated model is capable of providing equally as 
accurate Dielectric Function values as a non-correlated model. 
Correlation problems, however, are troublesome when arriving at 
the model, and when, for instance, attempting to apply it to 
alloy materials etc.) 

As regards application in the presently disclosed invention, 
it is important to realize that use of 983 Oscillator Structures 
can be difficult in conventional fitting precedures which require 
appropriate positioning of Oscillator Structures. The present 
invention overcomes this difficulty by simply defining segements 
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in an (e2) plot, and mechanically assigning an oscillator which 
has definite beginning and end points structure in each. 

Oscillator Structures are best described graphically. With 
5 this in mind it should be appreciated that the present invention 
will be better understood by reference to the Detailed 
Description Section, in coordination with the accompanying 
Drawings . 

10 
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BRIEF DESCRIPTION OF THE DRAWINaS 



Fig* 1 ^ows a plot of the imaginary component of an assumed 
Dielectric Function vs. Energy level, presenting with two Peak 
Regicms. 

Fig. 2 shows a basic Narrow Lor entzi an Oscillator Structure in 
which the tail regions are negligible. 

Fig. 3 shows the plot of Fig. 1, with two Narow Lorentzian 
Oscillator Structures applied thereto at the Peak Regions. 

Figs. 4a - 4d show various Oscillator Structures constructed by 
various application of Finite Order Polynomials and 
Discontinuities. 

Fig. 4e shows a Zero-Width Oscillator Structure, Cie. a Pole), 
constructed by application of coincident Lower, Center and Upper 
Energy levels and an Amplitude, or alternatively, a Discontinuity 
at a Center Energy level. 

Figs. 4f and 4g show Lorentz and Gaussian Oscillators. 

Figs. 4h, 4i , 4J and 4k respectively show lonicl, Ionic2, 
Harmonic, TOLO, and the like Oscillator Structures. 

Fig. 5 shows a plot of an actual Imaginary Component of a 
Dielectric Function vs. Energy level for a Gallium-Arsonide 
Sample System. 

Fig. 6 shows the plot of Fig. S with seven C7> Oscillator 
Structures applied thereto. 
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Fig. 7 demonstrates that the real part <el) of the dielectric 
function over an energy range of 1 to 6 eVy and demonstrates the 
presence of a Pole at about 6.5 ev- 

Figs. 8a - 8j demonstrate steps in the presently disclosed method 
of fitting Oscillator Structures to the imaginary part Ce2) of a 
sample Dielectric Function. 

Fig. 8k shows use of one Gaussian or Narrow Lorentz C6 or L), 
five Triangular <T> , and two modified triangular <T'> CT" ) 
shaped Oscillator Structures. 
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DETAILED DESCRIPTION 

It is announced at this point that Figs. 1 - 4e and 5 and 6 
are incorporated from Patent 5,796^983 for general insight^ Pigs, 
5 4f - 4k identify additional oscillator structures which can be 
applied in present invention^ and Figs. 7 - 8k demonstrate the 
presently disclosed invention methodology. 

As disclosed in Patent No. 5,796,983 and included herein for 

10 background and insight, there is shown in Fig. 1, a plot of data 
corresponding to the Imaginary Component (e2) of a presumed 
Dielectric Function vs. Energy, which plot is representative of 
that which might be obtained from investigation of a 
semiconductor material. Note the presence of two dependent 

15 variable "Peak" regions, located at "El" and "E2" on the 

independent variable "X" axis. Next, Fig. 2 shows the general 
shape of a Narrow Lorentzian or Gaussian Oscillator Structure 
centered about a point labeled "EC" on the independent variable 
"X" axis. (While not shown other Mathematical "Oscillator 

20 Structures" could also be utilized for the purposes of this 
demonstration). Note that the Narrow Lorentzian or Gaussian 
Oscillator Structure in Fig. 2 presents with a Maximum dependent 
variable Magnitude "A" and is shown with a "Full-Width-Half - 
Maximum" (FWHM) Broadening "B", (which Broadening can extend from 

25 Zero (0) to, theoretically, any finite value). Fig. 3 shows two 
such Narrow Lorentzian or Gaussian Oscillator Structures 
positioned on the (e2) vs. photon energy plot of Fig. 1. One 
said Narrow Lorentzian or Gaussian Oscillator Structure is shown 
with a Maximum Magnitude of "Al", a (FWHM) Broadening of "Bl", 

30 and it is centered about a point on the independent variable "X" 
axis labeled "El". The second Narow Lorentzian Oscillator or 
Gaussian Structure is shown with a Maximum Magnitude of "A2", a 
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(FWHM) Broadening of and It Is centered about a point on 

the independent variable "X" axis labeled "E2«. Note that the 
independent varaible "X" axis points "El" and "E2" in Fig. 1 and 
Fig. 2 plots correspond to one another. Fig. 3 also provides, 
5 in bold, a line representing the summation of the contributions 
from the two present Narrow Lorentzian or Gaussian Oscillator 
Structures, at each energy, which bold line very roughly 
corresponds to the (e2) vs. photon energy plot in Fig. 1. Note 
that the plot in Fig. 1 is, however, very different from the bold 

10 line plot in Fig. 3, particularly in the "Tail" regions beyond 
each "Peak", region. This demonstrates the shortcommings of an 
attempt to use Lorentzian and Gaussian etc. Oscillator Structures 
to model the Imaginary Component (e2) of a Dielectric Function 
vs. Photon Energy. The Three (3) degrees of freedom, (ie. 

15 Maximum Magnitude, Broadening and Center Energy), which define 
each Lorentzian or Gausian Oscillator Structure, are simply 
typically not sufficient to allow achieving a "good" Model fit in 
both Peak and Tail regions of most (e2} vs. Energy plots. 
However, note that sufficient data is contained in six (6) 

20 Lorentzian Oscillator Structure defining parameters to provide a 
rough approximation to the (e2) vs. Energy data. This compares . 
with hundreds, or even thousands of data points were the plot in 
Fig. 1 is determined and characterized as a table of data. Of 
course, additional Lorentzian Oscillator Structures could be 

25 added to the plot in Fig. 3, in an attempt to improve the 
approximation result, but this approach has not been found 
capable of providing results sufficient for use in Ellipsometer 
System settings, emphasis added. 

30 Note that Figs. 4f and 4g show Lorentz and Gaussian 

Oscillators in better detail. Shown are Central Energies (En), 
Amplitude (Amp) and Broadening (Br). It is to be appreciated 
that where the Broadening (Br) is small in the Lorentz Fig. 4f 
Oscillator, the Tail Regions quickly decay to an effective 



intersection with the horizointal axis, as demonstzated in Figs. 2 
and 3. Note, variations on the Lorentz and Gaussian Oscillators 
include lonicl, Ionic2, Harmonic, TOLO Structures as indicated in 
Pigs. 4h, 4i, 4j and 4k respectively. Said additional Oscillator 
Structures are described in the J. A Wool lam Co. WVASE Manual, 
which is incorporated by reference herein. Equations presented 
in the J. A. Woollam CO. WVASB Manual are included directly to 
provide insight to the Form of the Mathematical Equations which 
define them: 



Gausshm: 





Style 


EqoatiiM 


VHFmnaKata 


15 




-usftig 8S definod below. 




20 


GaujO(eV) 
GnLSCcm*^ 




Ampn ^Au((Bmentioiiless), 
Ent^E, (em''). BrI « Br, (em 
') 


25 


Cku.l(eV) 
GmjSlCMtf*} 


Br, Br, 


Em » &ref9. Bm-'Br^(eV) 

Ampn ^ Am (caf^), 

Em » Em tcm^)* Bm^^BK, (cm 

') 


30 


Gni2(eV) 
0811.7(011*') 


* Br, 


Ampn '=*Aa (dbtmdotdessX 
Em^En(elO.^°Br^(ey) 

Ampw^AmftSmea^n^s^ 
Emu ^ Em (em-'h Bm^Orm (att 



30 



Loiwrtz: 







FttParamctcfs 


LorjO(eV) 




Ampa ° 4, fiSawmimtlesaJ. 
Ent^E, (em-'). BrI - Br, (eaf') 




Lor.l (eV) 




Aapn '* AffeV). 

Em - E, (eV), Om'^^Br^ (eV) 

Bim "I^fem''). BrwBn, (an') 




|jw;2(eV) 




Bam^E,(eV),Bm^ar,(9V) 

Ampm^Am(itMnaoiies^. 
. Emi « E;,(em''). Em Br, (em') 





Harmonic: 



Style 


Eqnatioa 


FftFunmctars 


LorJO(eV) 




/lA^Mt <»/fa (tSnmasionhss^. 
Eni -^Em(an*h Bri « Br^feaf') 




ImA (cV) 
Um6(cin*') 




Sm » ^ ant » Ar. ^el9 




Lor^(eV) 
Lor.7(ciii'') 


AE/ 


A^m ^An (dbmndmUesa^^ 
Enn « i^re»9» Bam » Br. (eV) 

Amgm '^AmOSmemimUess). 
Em^Em (cm*). Om^Bir^ (eaf*) 


E/^E^^V4Br/^iBr^E 



31 



Ikinfr-f A fDn^c2r __. ^ n 


Style 




Fit parameters 


ioiilJO(eV) 
l(mi^(cai'') 




ecfcii " JB!jai (dSmensionles^ 
ewfii = ^ ((BSmemioabgaX 

tfcfcn » ^ifm (dbaenskudess), 

Eian » Er^ (cnf*h Bm^Br^fcm ) 


loii2.0(eV) 






£toit » (eaf'}. Bnt^Br„ (em ) 







1 \^M-X^. 


Fit Panundcn 


TOLOjO(cV) 
TOLO^(Gm'*) 




Elm o fteCe*?. ^ A<*« 

Ampa «^ frfimwiftwifcut. 

itteii - Btm (eat'). Bum = £U (coT) 


T0LO.i (cV) 
TOLO^tcm*') 




AmpH'^AmCy)- 

Ebm ^Etm<e^ ^ " i"**? 

Ebm "Etmfcm'). Eton "E^tcaTJ 
BtowBt^faa'). Aran -Ba.^ ) 


TOLO^CeV) 
T0L0.7(an*') 


, £^*-£'-iiWfi 


Bom''EtmftV),Etm''Eum(eV) 

itinpa Am/dlmeaaoakss). 

Ekm^Etmfeaf'i.Eum-^E^iem'j 

BtowBi^ibarj.Sum'Bmfar) 



32 



It Is to be appreciated that for application in the 
presently disclosed invention, the functional requirement is that 
an Oscillator must provide substantially definite starting and 
ending intersection points with the horizontal axis so that an 
5 end point of one Oscillator can be vertically aligned with peak 
point of an adjacent Oscillator. The 983 Patent Oscillator 
Structure which is reviewed herein in detail, is preferred as it 
can be constructed to meet said functional reguirment as 
demonstrated well in Figs 4a - 4c. 

10 

The forgoing then, demonstrates the concept of utilizing a 
summation of contributions from a plurality of Mathematical 
Structures to attempt to Model a plot of, for instance, (e2) vs. 
Energy data, in a manner which requires determination of a 

15 reduced number of parameter values, (as compared to an approach 
utilizing, for instance, tabulated data), and the foregoing also 
identifies the shortcommings of the demonstrated approach. The 
preferred realization of the present invention serves to overcome 
the identified shortcommings by utilizing novel Oscillator 

20 Structures which. provide more available degrees of freedom than 
provided by, for instance, the Lorentzian Oscillator Structure. 
This is accomplished by utilization of defined novel Oscillator 
Structure(s) , which allow shaping the "Peak" and "Tail" regions 
thereof, (which "Tail" Regions project to the right and left of a 

25 "Peak" Region Center Point thereof), in a very non-correlated 

manner. The use of a "Broadening Factor" applied over the full 
range of an Oscillator Structure, conceptually similar to that 
described infta with respect to the Lorentzian Oscillator 
Structure, is retained, but application thereof to a present 

30 invention Oscillator Structure is accomplished via Mathematical 
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Convolution Integrations applied simultaneously, but 
independently, to various regions of an Oscillator Structure. 
This is required, because as described directly herein, different 
independent variable regions of an Oscillator Structure are 
defined, for instance, by contributions from different Finite 
Order Polynomials and/or Discontinuities. 

Turning now to Fig. 4a, there is shown a general 
representation of a novel Oscillator Structure disclosed in 
Patent Mo. 5,796,983. Note that Four (4) Finite Order Polynomial 
defined regions, and one Finite Magnitude Discontinuity are 
shown. The Four (4) defining Polynomials are identified as PI, 
F2, F3 and F4, with Fl beginning at coordinates (0,EL), and 
merging with Polynomial F2 at coordinate (ALM,ELM}. Polynomial 
F2 is terminated at coordinate (AL,EC) by the lower end of the 
LDISC region of the Discontinuity, (shown as two dimensionless 
regions, labeled UDISC and LDISC, above and below coordinates 
(A, EC)) respectively. (Note that the magnitudes of UDISC and 
LDISC are typically selected to be of equal finite Dimensionless 
Magnitudes). The upper end of the UDISC region ends at 
coordinates (AU,EC), whereat Polynomial F3 begins. Polynomial F3 
extends to a merger with Polynomial F4 at coordinates (AUM,EUM), 
and Polynomial F4 intersects the independent variable "X" axis at 
coord i na te ( 0 , EU ) . 

As required Oscillator Structure defining Mathematical Model 
Coefficients for the above described Oscillator Structure (s ) , the 
following, underlined . Nine (9) Coefficients were chosen in the 
983 Patent disclosure for evaluation and application in the 
practice of the Method of use described therein: 
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AU = A*(1+UDISC); 

AL = A*(1-LDISC); (Note UDISC and LDISC are of equal 

magnitude) 

AUM = AU* UainD ; 

ALM = AL* Lamp ; 

EUM = EC + (EU-EC)* Upos ; 

ELM = EC - ( EC-EL ) * Lpos . 

In addition, a Broadening Factor, Sioma (O" ), assumed to 
be, but not necessarily, Gaussian, is also evaluated to determine 
a defined Oscillator Structure. (Note, Gaussian Broadening is 
preferred because convenient, and novely applied. Mathematical 
aides are available when such is the case, as described elsevhere 
in this Disclosure. It is to be understood that other Broadening 
criteria can be utilized and remain within the scope of the 
present invention.) 

The "Order" of the various Finite Order Polynomials must 
also be identified, and such is accomplished with the use of two 
(2) additional parameters as follows: 

For EL<E<ELM: 

Fl = ALM*( (1-L2d)*yl + L2d*yl^2) 
For ELM<E<EC 

F2 = (ALM + (l-ALM)*(l-cl-dl)*y2 + cl*y2"2 + dl*y2"4)) 
For BC<E<BUM 

F3 = (AUM + (l-AUM)*(l-cu-du)*y4 + cu*y4'^2 + du*y4'"4)) 
For EUM<E<EU 

F4 = AUM*( (1-U2d)*v3 + U2d*y3^2) 

where: 

Cl = L2d* (ALM/ ( 1-ALM) )*( (EC-ELM) /(ELM-EL ) )'^ 2 
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dl = {(1/(1-ALM))*(1-((EC-ELM)/(ELM-EL))*ALM*( EC-ED* 
( ( L2d/ ( ELM-EL) + (1/(EC-BLM)) 

cu = U2d*(AUM/(l-AUM) )*( (EC-EUM)/(EUM-EU) )"2 

5 

du = ( (l/d-AUM) )*(!-( (EC-EUM)/(EUM-EU) )*AUM*(EC-EU)* 
( ( U2d/ (EUM-EU) + (1/{EC-EUM)) 

and where: 

10 

yl = (E-EL)/( (ELM-EL); 
y2 = (E-ELM)/ (EC-ELM); 
y3 = (EU-E)/(EU-EUM); and 
y4 = (EUM-E)/(EUM-EC); 

15 

(Note that the yl and y2 as used here are not the same as the 
yl and y2 as defined with respect to EQ, 3 in the Disclosure 
of the Invention Section herein. The yl and y2 in the 
disclosure of the Invention Section are related to a Gaussian 
20 Broadening Factor ) • 

are normalized energy terms. The significance of the use of 
"normalized" energy terms is particularly critical in application 
where Gaussian Broadening of Oscillator Structures is utilized. 

25 Briefly, use of normalized energy terms allows Convolution 
integration of a product of a Finite Order Polynomial and a 
Gaussiam Broadening Factor to provide one-dimensional, 
normalized, independent variable based Look-up Tables to allow 
convenient evaluation of corresponding dependent variables. Said 

30 significance is described more fully in the Disclosure of the 
Invention Section herein. 

Fig. 4b shows another Oscillator Structure formed from the 
same basic components as that in Fig. 4a, (eg. Finite Order 
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Polynomials and Finite Magnitude Discontinuities), but in which 
Polynomial F2 and the upper end of UOISC meet at coordinates 
(AU,EC) and Polynomial F3 begins at coordinates (AL,EC) and ends 
at coordinates (0,EU). Note that no Polynomial F4 is present. 
5 This is demonstrative of the fact that some or all of the Four 
Polynomials shown in Fig. 4a can be absent while still defining 
an Oscillator Structure. In that light note that Fig. 4e shows 
the case where no F1-F4 Polynomials are present, with the result 
being a '•Zero-Width-Finite-Magnitude Oscillator Structure", or 

10 "Pole", of Zero-Width, (ie. (EL), (EC) and (EU) are coincident), 
but with a Finite Magnitude (AU) above the independent variable 
"X" coordinate (EC). (Note that Fig, 4e can alternatively be 
considered to show a Discontinuity of Finite Magnitude (AU) above 
the dependent variable zero (0) level, present at photon energy 

15 (EC), said Finite Magnitude Zero-Width Oscillator and Finite 
Magnitude Discontinuity interpretations being equivalent 
mathematically). Fig. 4c shows a case wherein Polynomials Fl, 
F2, F3 and F4 are present, but no Discontinuity DISC, (ie. UDISC 
and LDISC, which are of equal Finite Magnitude), is present. 

20 Polynomials F2 and F3 merge at coordinates (A, EC), which 

correspond to a "Peak" in this Figure. Fig. 4d shows a case in 
which a Discontinuity is present at "X" axis coordinate (EC), and 
in which Polynomials F3 and F4 are present, but Polynomials Fl 
and F2 are absent. It is to be understood that an analogous case 

25 wherein a Finite Magnitude Discontinuity and Finite Order 

Polynomials Fl and F2 are present to the left thereof, but where 
Finite Order Polynomials F3 and F4 are absent can also be 
constructued and such consitutes an Oscilator Structure. It 
should also be appreciated that a single Finite Order Polynomial 

30 with comprising Terms thereof of opposite signs can serve to 
provide a peaked Dependent Variable vs. Independent Variable 
Plot, which begins at an dependent Variable value of zero (0.0), 
and returns to a dependent Variable value of zero (0.0) at some 
Dependent Variable value. For instance the Finite Order 
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Polynomial (Y = 3*X - begins at a "Y" value of zero (0.0) at 

"X" = (0.0) and returns to a "Y" = (0.0) at "X" (3.0). Such a 
Single Finite Order Polynomial can serve as an additional 
Oscillator Structure^ (over a range of (0.0) to (3.0)), although 
coordination with other Oscillator Structures and application of 
Convolution Integral effected Broadening thereto, presents 
difficulties. 

It is to be noted that where two Polynomials merge, (eg. FX 
and F2, F3 and F4, or, possibly, when present without a Finite 
Magnitude Zero-Width Discontinuity therebetween, F2 and F3), they 
do so with continuous zero, first. and preferably second (2) 
derivatives at the point of their merger. This is handled by the 
Mathematical Model Coefficient evaluation programming involved in 
implementing the invention. 

It is also noted that the straight line components of the 
Oscillator Structure shown in Figs. 4a - 4e can be combined to 
form Triangular (T) Oscillator Structures as demonstrated in 
Figs. 8a and 8b, and shown applied in Figs 8c - 8k. Said 
Triangular (T) shaped Oscillator Structures are a preferred, 
though not exclusive, structure in practice of the present 
invention methodology. 

It should then be appreciated that the Oscillator Structures 
just described be constructed from one or more Finite Order 
Polynomial (s) and Finite Magnitude Discontinuities in various 
sequences, and that in the presently preferred embodiment, a 
typical, non-limiting Oscillator Structure can include from zero 
(0) to Four (4), Finite Order Polynomials, each present Finite 
Order Polynomial including terms of any Finite Order, (eg. zero 
and/or first and/or second and/or third etc.), and being 
positioned to the right or to the left of a Center Point; and 
that Oscillator Structures can include Finite Magnitude 
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Discontinuities therein, present at Center Points thereof. It is 
also to be understood that acceptable such Oscillators include 
those of Zero (0) width, (ie. Poles) in which no Polynomials are 
present, said Poles being mathematically equivalent to a single 
Zero-Width Finite Magnitude Discontinuity, beginning at a 
dependent variable value of zero and positioned at a specified 
independent location. 

In conventional useage, various constructed Oscillator 
Structures are positioned, by a user, at appropriate locations in 
regions along an independent variable "X" axis, which locations 
correspond to "Peak" locations regions in a plot of, for 
instance, the dependent variable Imaginary Component of a 
Dielectric Function (e2) vs. Photon Energy, and via Mathematical 
Convolution Integration, appropriate constructed Oscillator 
Structures are subjected to Gaussian Broadening. 

Turning now to Fig. 5, there is shown a plot of an actual 
Gallium-Arsonide Sample System Imaginary Dielectric Function (e2) 
Component vs. Energy. Fig. 6 shows the plot of Fig. 5 
approximated by the Summation of the contributions from seven (7) 
Oscillator Structures as disclosed in Fugs. 4a - 4e, (one of 
which, (ie. Oscillator Structure #1 in following Table 1), is not 
readily visible, but is present at the left side of the plot, 
with an "X" coordinate (EC) of 1.785 ev). Fig. 6 demonstrates 
the conventional approach to using said Oscillator Structures to 
fit a Complex Curve. Note that each of the Seven (7) Oscillator 
Structures is applicable over only portions of the full Energy 
range in the Fig. 6 plot. The Finite Order Polynomial Coefficient 
and Finite Magnitude Discontinuity definign Mathematical Model 
Coefficients, (see infra for definition thereof), which serve to 
define the seven (7) Oscillator Structures utilized in 
constructing the Mathematical Model represented in Fig. 6 are 
provided in Table 1. 
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A Zero-Width Finite Magnitude Pole with a Magnitude of 50.75, 
present at 7.195 ev, (beyond the range of Energy levels shown in 
Fig. 6), tabulated in the seventh row of Table 1, was assumed as 
present to account for all effects caused by unknown data to the 
right of 6.5 ev. The value of the Magnitude thereof being 
determined by the parameter evaluation routine which also 
provided values for the various Finite Order Polynomial 
Coefficient and Finite Magnitude . Discontinuity Magnitudes 
defining Mathematical Model Coefficients. (Note, a 
mathematically equivalent result could be effected by use of a 
Finite Magnitude Discontinuity of 50.75, beginning at at 
coordinates (0.0, 7.195), rather than by the use of a Zero-Width 
Oscilator Structure, (ie. a "Pole"), as shown in Table 1). 



Comparison of Figs. 5 and 6 show that the Mathematical Model 
arrived at by practice of the Method of Use of the Oscillator 
Structures as applied to the plotted data in Fig. 5, provides a 
very close replication of said data, as shown in Fig 6. With 
that in mind,, it is noted that only Oscillator Structure 2, (see 
Table 1), provides second order Polynomials, (ie. L2d is 1), and 
that Oscillator Structures #0 and #1 have no associated Fl and F2 
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Polynomials present. It should be appreciated that the 
Oscillator Structure Mathematical Model Coefficients identified 
as determined in Table 1 and Fig. 4a, are arbitrary, and that any 
Coefficients which allow evaluating present Finite Order 
5 Polynomial Coefficients and Finite Discontinuity Magnitudes, 
consistant with application of Convolution Integral effected 
Oscillator Structure Broadening, could be selected as appropriate 
for determination determination and use in calculation of 
dependent variable values given independent variable values, 
10 which in the above provided example were (e2) vs. Photon Energy. 

It is noted at this point that Fig. 6 demonstrates typical 
conventional application of the 983 Patent Oscillator Structures 
to model Dielectric Functions. Such conventional practice 

15 strives to obain good fit with a minumum of Oscillator 

Structures, each of which has an effect over a fairly large 
energy range and can be difficult to apply, even for experienced 
practitioners. In comparison, as will be discussed with respect 
to Figs. 7 - 8e, the presently disclosed invention makes use of a 

20 multiplicity of short-energy-range effect oscillator structures 
to the end that the imaginary part of a Dielectric Function is 
effectively modeled as a sequence of line segements between 
adjacent Oscillator Structure Peaks, the shape of said line 
segements each being determined as a sum of the Oscillator 

25 structure components present between said adjacent peaks. 

Turning now to Fig. 7, which introduces the approach of the 
presently disclosed invention, it should be noted that it 
demonstrates that the real part (el) of the dielectric function 
30 over an energy range of 1 to 6 eV, and demonstrates the presence 
of a Pole at about 6.5 ev. An insert shows that a model for the 
sample for which the can comprise a substrate (SUB) and a Thin 
Film (TF). The presently disclosed invention teaches that the 
Thin Film (TF) can be mathematically modeled by a Thickness and. 
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for instance^ a Cauchy Dispersion Function and parameters 
therein^ along with Loation and Amplitude for the Pole (P) 
evaluated by a global regression of said model onto the real part 
(el) of the dielectric function data. (See Fig. 4e for insight 
5 as to how a Pole can be formed by setting Fl, F2, F3 and F4 to 
zero in Fig. 4a). A first fit step in the presently disclosed 
invention is then to determine Thin Film (TF) Thickness and 
parameter values which define the Pole (P) location and amlitude^ 
(which pole accounts for effects at energies higher than a 
10 defined energy range of interest). 

Figs. 8a ~ 8f demonstrate the steps of the presently 
disclosed invention which involve a method of applying a 
plurality of applying oscillator structures which are formed by 

15 variously assigning definition to Fl, F2, F3 and F4 in Fig. 4c, 
(ie. Fl and F4 are assigned to be straight lines, and F2 and F3 
set to zero lenght). Note that a similar oscillator structure 
results from setting F2 and F3 to be straight lines and Fl and F4 
to be of zero length, (again see Fig. 4c for the locations of Fl, 

20 F2, F3 and F4). Other combinations are possbile to provide the 
same functional result. An example of the resulting oscillator 
structure can be observed in Fig. 8b, where is shown a triangular 
shaped sturcture in dashed lines spanning energy 3 - 5.5 eV with 
a peak at 4.5 eV. Note that Fig. 8b also shows that additional 

25 oscillator structures are located to the right and left of the 
oscillator centered at 4.5 eV, and that a common theme is that 
the disclosed oscillator structures begin and/or end where an 
adjacent oscillator structure peaks. The' Triangular Structures 
formed by the straight lines, (eg. Fl and F4), are each 

30 identified by (T). Figs. 8a and 8b show simple demonstrative 

imaginary dielectric part (e2) plots which are fit, respectively, 
by two and three. Triangular SHaped Oscillator Structures (T) . 

Figs. 8c and 8d demonstrate application of the method of the 
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presently disclosed invention. Fig, 8c shows a first Oscillator 
Structure (OSl) applied to fit a left-most portion of a shown 
imaginary part of a dielectric part (e2), and Fig. 8d shows first 
and second Oscillator Structures (OSl) applied to fit a bit more 
of the left-most portion of said shown imaginary dielectric 
function part (e2). (Note, the procedure could have been started 
at the right side of the plot and remain within the scope of the 
disclosed invention). In practice the first Oscillator Structure 
(OSl) is fit to the data, typically via a point-by-point square 
error reducing procedure. This is followed by adding the Second 
Oscillator Structure (0S2) as shown in Fig. 8b, and again 
performing a fitting procedure which evaluates parameters in said 
second oscillator structure (0S2) and simultaneously re-evaluates 
parameters in the first Oscillator Structure (OSl). Fig. 8e 
demonstrates a result where the point-by-point procedure 
indicated by Figs. 8a and 8b is continued over the entire energy 
range. Note that Fig. 8e assumes that equal energy (ie. 
wavelength), range segement steps are used for all oscillator 
structures. It can be observed that this leads to error in peak 
regions. A modified approach is demonstrated in Fig. 8f, and 
provides for using variable length wavelength range segements/ 
specifically note the use of smaller energy steps in peak 
regions. Again, as Oscillator Sturcture are added a fitting 
procedure is again performed which evaluates parameters in an 
added oscillator structure and simultaneously re-evaluates 
parameters in the previously evaluated Oscillator Structure (s ) . 

With all Oscillator Structures evaluated by the 
point-by-point approach just described, using parameter values as 
starting values, a global fit can be performed. It is noted that 
typically the evergy range over which an Oscillator Structure 
applies is fixed during point-by-point procedures, but in such a 
global fitting, the base width as well as amplitude of the 
Oscillator Structures can be allowed to float. Further, while 
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the point-by-polnt fits are onto data representing the imaginary 
part of the Dielectric Function, the Global fit can be onto data 
representing both real and imaginary parts of thereof. 

Figs. 8g - 8j demonstrate that the approach to assigning 
oscillator structures can involve beginning at a central location 
and specifying some number of segements. Fig. 8g speficies one 
segement, while Figs. 8h and 8i specified two and four equal 
length segements respectively. Fig. 8j shows five specified 
segements, wherein the last two are shorter than the first three. 
Note that where two segements are specified, three oscillator 
structures are present. Fig. 8i shows that where four segements 
are present, seven oscillator structures are present and Fig. 8j 
shows that where five segements are present nine oscillators are 
present. This occurs based on the method of positioning the 
oscillator structures such that sequentially, each successive 
oscillator structure begins where the just prior oscillator 
structure peaks. Typically global fitting will be applied in 
such a scenario. 

Fig. 8k shows use of one Gaussian or Narrow Lorentz (G or 
L), five Triangular (T), and two Modified Triangular (T') (T' • ) 
shaped Oscillator Structures. Note that Modified Triangular 
Shaped Oscillator (T*) includes Discontinuies (UDISC) (LDISC), 
and that Modified Triangular Shaped Oscillator (T* • ) is 
constructed from a polynomial (Fl) and a straight line segement 
(F2). Fig. 8 is included to demonstrate that while triangular 
(T) shaped oscillators are preferred, the present invention 
methodology allows for use of different functional Oscillator 
Structures in modeling a dilectric function. A requirement of 
all Oscillator Structures is, however, that definite start and 
end point intersections with the horizontal axis must be 
identifiable so that generally the start and end points of one 
oscillator can be positioned directly beneath the peak of the 
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Oscillator Structures directly adjacent thereto. It is pointed 
out that the start and end of the first and last Oscillator 
Structure, of course, have no adjacent Oscillator Structure peak 
positioned thereabove. 

It is emphasised that while, for simplicity. Triangular 
shaped Oscillator Structures are shown in Figs. 8a - 8j, any 
functional shape Oscillator Structures such as shown in Figs. 4c 
and 4f - 4j can be applied, and it is possible to apply more than 
one type, (eg, mixed Traingular and/or Modified Triangular and 
Gaussian and/or Narrow Lorentz and/or Harmonic and/or Triangular 
and/or TOLO etc., see Fig. 4k for an example). The novel 
requirement is that the segmental modeling of imaginary parts of 
dielectric functions is with Oscillator Structures that have 
substantially starting and ending points and which are 
Kramers-Kronig (K-K) consistent be achieved. 

It is be understood that while Dielectric Functions are 
recited in the Claims, Refractive Index and Extinction 
Coefficient Optical Constants are mathematically equivalent and 
in that light are within the scope thereof. 

It is also to be understood that the terminology "Pole" can 
include functional equivalents to zero width structures. For 
instance, very narrow oscillators with finite amplitudes are 
within the scope of "Pole". 

It is also to be understood that providing experimentally 
obtained data for real and imaginary parts of a dielectric 
function vs. wavelength for a sample comprising a transparent 
thin film on a substrate, over a determined range of wavelengths 
can involve providing equivalent ellipsometric PSI and 
ellipsometr ic DELTA data. 
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It is to be understood that Dielectric Function vs. Photon 
Energy data used in the practice of the Method of the present 
invention can be obtained by Ellipsometry, and by any other 
applicable technique, such as Transmission and Reflectance. 

Having hereby disclosed the subject matter of the present 
invention, it should be obvious that many modifications, 
substitutions, and variations of the present invention are 
possible in light of the teachings. It is to be understood that 
the invention can be practiced other than as specifically 
described, and should be limited in breadth and scope only by the 
Claims. 
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